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ANALYSIS OF DEVELOPING LAMINAR FLOW AND HEAT TRANSFER 


IN A TUBE FOR A GAS WITH COOLING 
by Alden F, Presler 
Lewis Research Center 

SUMMARY 

The compressible boundary layer equations are solved numerically for subsonic 
flow of a monatomic gas in a circular tube with uniform cooling flux on the walls, 
hiitial fluid conditions are uniform, with velocity and temperature profiles developing 
downstream of the tube entrance. 

Both moderate and high cooling flux conditions are investigated, including their 
effect on profiles, static pressure development, Nusselt numbers, and friction param- 
eters based both on local shear stress and pressure drop. 

Heat transfer and friction factor values for flows at moderate cooling approach 
those found in the constant property model. High cooling rates, however, distort the 
velocity profiles to such an extent (due to adverse pressure gradients) that Nusselt 
numbers and friction parameters deviate greatly from the standard constant property 
model. 

Downstream of the entrance region, a very simplified calculation of the momentum 
pressure drop with Rayleigh flow assumptions gives values within 25 percent of the more 
exact boundary layer results in the case of the higher cooling flux. 


INTRODUCTION 

Reference 1 provided one of the earliest analyses of laminar flow with cooling which 
allowed temperature variation of the physical properties. The analysis of reference 1 
was essentially fully developed in that the velocity and temperature profiles were inde- 
pendent of the axial coordinate. Heat transfer and friction parameters were obtained for 
both cooling and heating cases. 

Reference 2 reported the numerical analysis of the compressible boundary layer 
equations for low subsonic gas flow with developing velocity and temperature profiles 



and for both high heat flux input and essentially adiabatic conditions. Rather large devia- 
tions in fluid friction and heat transfer parameters compared with constant property re- 
sults (ref. 5, chs. 9 and 13) were shown for cases of high heat flux. 

The analytical portion of the present report has its basis in reference 2, and, in 
extending that into the problems of cooling heat transfer, it details the flow history of 
static pressure development and its momentum and frictional components. 

This report contains a finite difference numerical analysis of the subsonic laminar 
flow of a compressible gas in a tube with uniform flux cooling heat transfer. Monatomic 
gas properties are used. The boundary layer equations are integrated across the tube 
radius without the use of the inviscid core assumption. The flow entering the tube is 
uniform in temperature and velocity, and these profiles develop with the flow through the 
tube. Effects of the level of cooling flux on local static pressure and velocity and tem- 
perature profiles are obtained, and from these basic variables are derived the variations 
in heat transfer and friction parameters. In this analysis all physical properties except 
specific heat have temperature dependency. 


ANALYSIS 


The steady-state momentum and energy boundary layer equations and the continuity 
equation, in axisymmetric cylindrical coordinates, are, respectively (rt... o) 


9z 9r dz r 9r \ 9r/ 


( 1 ) 


pc 


P 




dz r 9r \ 9r/ \9r/ 


( 2 ) 


and 


^ + (3) 

9z r 9r 

(Symbols are defined in the appendix. ) The perfect gas law is 

p = p^t (4) 

The boundary and initial conditions for the system are 
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Initial (z = 0, 0 < r < r^): 


V = 0 
P=Pi 


Centerline (z > 0, r = 0): 


u = Uq(z) 

V = 0 

t = tQ(z) 

3u _ 0t _ Q 
dr dr 


Wall (z > 0, r = r^): 


u = V = 0 


+ 





= Constant 


The initial conditions are for uniform temj^erature and velocity profiles. The cen- 
terline conditions are results of the axial symmetry assumption. The wall conditions 
reflect the nonporous wall statement and the thermodynamic convention that heat added 
to the system is positive. Figure 1 is a model of this sytem., including initial and 
boundary conditions. 

The system of equations (1) to (4) and the boundary conditions may be nondimen- 
sionalized by means of the following terms; 
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The nondimensional system with boundary conditions is 


A (rm'^ 

V az dR dz R aR V eR 


p'CpIu ^ + V 

n az 3R 


axA (y - _ dp+ 4 1 


u 


A (kr 


+\ 4(y 


Z Pr. R aR 


aR 


AA + 1J_ (p'VR) = 0 
az R 3R 


1 - yMfP'^ 

p’ = 


Initial conditions (Z = 0, 0 ^ R ^ 1): 


1 - q+T+ 


U = 1 
V = 0 


P"^ = 0 


T+ = 0 


Centerline conditions (Z > 0, R = 0): 



u 


= u. 
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Wall conditions (Z > 0, R = 1): 


U = V = 0 


/3T'*~\ ^ J_ 

\0R/^ 

K(Z, 1) = 
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The radial velocity term for both energy and momentum equations is obtained in 
terms of the axial velocity from the continuity equation 


P’ 




dR 
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where 
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The axial pressure gradient is obtained from the momentum equation at the wall 
(R = 1, U = V = 0): 
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The thermal and transport properties are evaluated by the special relations 


M’ = K = 


ir 


(1 - q'^T 


+^ 0,68 





Pri = 


2 
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and 





(16) 


These correspond closely to the properties of helium gas (ref. 2), 

The radial derivatives of the transport properties are, for the preceding special 
relations, 


dll' ^ -0,68 q+ /aT~^\ 

0R, 0R \ / 


(17) 


which indicates that for vanishing heat flux we can assume constant property flow. 

Numerical solution of the boundary layer system is accomplished by replacing the 
axial and radial derivatives by the finite difference approximations: 

3U_ au _ ^j + l,k ~ ^j,k (- 

az az AZ 


^ ^j,k+l ^j,k-l 

aR AR 2AR 


(19) 
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( 20 ) 
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where j, k are the axial and radial indices, respectively, of the finite difference mesh. 


Friction and Heat Transfer Parameters 


The local pressure and pressure gradient are calculated at each axial position, and 
the temperature and velocity profiles also as a function of radius. With these basic 
quantities, certain useful integral quantities can then be obtained. 


= 2 f RUdR 
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In terms of the dimensionless temperature ratio 0 = (t - t^)/ (t^ - t^), the bulk Nusselt 
number is 


Nu, = • — 


( 29 ) 


Using the power relation for thermal conductivity, 


Nu^ = 2/^) ■ 

\tb 


, 0.68 


(30) 


The friction factor based on the static pressure drop from the entrance to the axial 
distance z/D is defined as 


P,b 




p) 


z \ 2 


(31) 


and in the dimensionless variables as 


p,b 


ip+ 




(32) 


or 


4f = 

p, b b 


2P~^ 


(33) 


The friction factor based on the local shear stress is defined as 


2t. 


T,b 


w 


^b^b 


(34) 
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and the shear stress parameter of reference 1 is 



Then, from equation (34) and (35) we obtain 


4f j^Re^ = 


( 35 ) 


(36) 


= 166 j, 


(37) 


The shear stress parameter is also written in dimensionless terms as 



(38) 


where 


^w 


T r 
w w 


MiUi 


(39) 


The numerical program provides values of P'*', Uj^, f^ j^, and 6|^ as functions 
of Z"^ and cooling flux q'*'. 


Momentum and Friction Pressure Drops 


The decomposition of static pressure drop into friction and momentum components 
is readily followed through the impulse function of compressible gas dynamics (ref. 4, 
chs. 5, 8). The impulse function per unit area is 

I = p + p^ul (40a) 

= p + Guj^ (40b) 
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where G is the mass flux. From the one -dimensional momentum equation 


dl - dp + G du^ 


= -4r. 


w 



( 41 ) 


where is the local wall shear stress. Therefore, on integration of equation (41) 

AI = Ap + G Auj^ (42a) 

= -47w(5) 

where L is the distance between measuring stations and is the average shear stress. 

Taking the lower limit of integration as the initial position or entrance to the tube, 
equation (42a) is 


AI = I - L 


(43a) 


AI = (p - p.) + G(Ujj - Uj) 


= -4? (h 

'"Vd 


(43b) 


where u. is the bulk velocity at the axial position Z units from the entrance. 

^ 2 
Dividing both sides of equations (43) by results in 


AI 


4t, 


w 


PfU, 


2 2 
p.U. 
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-) 

D/ 


(44a) 


AI Pj-P , -^j) 


PfUi 


PfUj 


p.u2 

^11 


(44b) 


Since mass flux G = Pj^Uj^ = ^i^i» Ihen equation (44) is written as 
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( 45 ) 
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w 
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Pi - p , 

Pi“i Pi“i 



= -P+ + (Ub - 1) 


Equation (45) states that the overall static pressure drop is the sum of the drop due 
to friction and the drop due to momentum change. Thus, 


Aprp = Appj^ + Ap 


•M 


(46) 


Divide both sides of equation (46) by PjU^, 


AP-P Appj^ Ap 


M 


2 2 2 

p.u. p.u. p.u. 

^11 ^11 '^11 


(47) 


the left side of equation (47) is identifiable as the dimensionless pressure P'*'. Hence, 
equation (47) is also written as 


where 




4t, 
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FR 


PiUi 



(48) 


(49) 


and 


Pm = (Ub - 1) (50) 

The numerical analysis provided values of P'*’ and 11^ as functions of z/D for 

each q'*' and M, case. Equations (48) and (50) are used to compute the frictional Pi,rj 
1 r K 

over the axial distance z/D for the one -dimensional compressible analysis. 
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Rayleigh Flow 


Frictionless one -dimensional compressible tube flow with heat transfer (Rayleigh 
Flow) can be readily calculated for the uniform heat flux case. This simplified treat- 
ment provides a convenient way of comparing the idealized momentum change with the 
actual momentum change in a cooling situation. 

The very complete and readable source for one-dimensional compressible flow 
theory in reference 4 is the background for this section and is taken from chapters 9 
and 20 of that reference. 

The energy balance for the tube with uniform heat flux and for the gas with constant 
specific heat is 


T — Q 
^2 


+ T, 


CpW 


(51) 


where positions 1 and 2 are separated by tube distance z. Let us take position 1 as the 
initial point of flow and heating. The preceding equation is nondimensionalized with the 
initial temperature L, 


that is, 


Again, 




C t.W 

p 1 



(52) 



(53) 



(54) 


which, for fixed Mj, tj, and q'*', is a linear equation in the reciprocal Graetz variable. 
The temperature ratio in equation (54) is further nondimensionalized with the critical 
stagnation temperature T* so that 
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( 55 ) 


X 

T* ^ 8q~^ • Gz ^ j 

T. ■ t! 


where for any Mach number 


■pt 


2(y + 1)M^ 1 + ILlA 
2 

[l + 


(56) 


Again, for fixed M-, t., and and for given values of Gz we calculate from 

J. X ^ 

equation (55) a succession of values of T/T . The Mach number M corresponding to 
this temperature ratio can be calculated from equation (56), or, more conveniently, from 
charts of T/T against M in reference 4. 

The Mach number variation with heating (cooling) length is now the prime variable 
in calculating the static pressure development in the tube. Since impulse function is 
constant for frictionless flows. 


P2 + 




(57a) 


or 


P2(l + rM2) = Pi(l+yM2) 

Then, 


P2 _ 1 + yM^ 

Pi 1 + yMg 


(57b) 


(57c) 


where Mach number is that of the local bulk velocity and sound velocity. If we identify 
Pj and Mj^ as belonging to the tube inlet. 
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Pj 

pA 


?2 

Pi 


1 


1 + yM^\ 

1 H- rmlj 


A - A 

Mj (l + yM|) 


(58) 


A simplified treatment is possible for flows with low subsonic Mach numbers (say 

2 

M < 0, 1). For this case equation (56) can be written as 


— = 2(y + 1)M^ 


and therefore, equation (54) as 




Also, equation (58) is approximately 



M? M? 


Combining equations (54a) and (58a) gives us 



(56a) 


(54a) 


(58a) 


(59) 


Equation (59) states that the dimensionless pressure drop for a frictionless gas flow, 
which is the same as the momentum pressure drop of equation (50), is directly propor- 
tional to the dimensionless flux and the dimensionless distance down the tube (inverse 

2 

Graetz number), provided that M < 0. 1. 
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Equation (59) predicts an adverse pressure gradient in all sections of the tube for 
cooling in the absence of friction. 


RESULTS AND DISCUSSION 

The system of equations (5) to (11) were solved numerically for radial variations of 
U and T and for the axial variation of U, T, and These dependent variables 

were used in further calculation of the heat transfer (Nusselt) parameter and of the 
various and useful friction parameters. Calculations were carried out for combinations 
of values of the initial Mach number, = 0. 004, 0. 06, and 0. 30, and of the heat flux 
parameter, q'*’ = -0. 135, -1.35. 


Velocity Profiles 

Radial profiles of the axial velocity components for several axial positions, initial 
Mach numbers, and flux parameters, are shown in figure 2. 

Development of the velocity profiles for the low heat flux cases (figs. 2(a) and (b)) 
follows closely that for laminar incompressible flow (ref. 5 - fig. 9.3, p. 231), becom- 
ing essentially a Poiseuille profile for distances of Z'*' > 0. 10. 

For the high cooling rate cases, the profiles for Z'*’ = 0. 1 reflect a cold dense gas 
along the wall which is decreasing in velocity under both adverse pressure gradient and 
density change, thus "thickening the boundary layer" and forcing more of the flow 
toward the centerline. The resulting centerline velocity excess in figures 2(c), (d), and 
(e) is opposite in behavior to the velocity profiles for gas cooling found in references 1 
and 6 (fig. 9. 16 therein). In reference 1 the radial velocity was omitted from the 
momentum equation, and the axial gradient 3u/9z was assumed to be small. These 
restrictions imposed a fully developed condition on the flow, under which the profile 
shape was determined by the temperature behavior of thermal conductivity and viscosity. 
Since gas viscosity decreased with cooling, shear stress decreased near the wall and 
velocity increased in that region thus decreasing the centerline velocity needed to main- 
tain a given bulk velocity. 


Temperature Profiles 

Figures 3(a) and (b) for low heat flux conditions show that the temperature profile 9 
development in the axial direction is like that of the velocity profile. Values of 6 for 
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= 0. 1 compare closely with fully developed values for t^/t^^ = 1 from reference 1 
in which the flattening of the profiles is due to increased heat conduction in the wall 
region because of the lowering of gaseous thermal conductivity with decreasing wall tem- 
perature. In the present analysis for higher cooling flux the thickening boundary layer 
also convects cold fluid into the fas moving center stream. This augments the cooling 
because of the lower thermal conductivity near the wall, the effect of which on tempera- 
ture profiles for = 0. 1 is seen in figures 3(c) to (e). These profiles are flatter than 
the fully developed case (fig. 3 of ref. 1). 


Nusselt Numbers 

Development of bulk Nusselt number with axial distance varies with the level of cool- 
ing heat flux. For initial Mach numbers of 0. 004 and 0. 06, figure 4(a) compares the 
values of Nu^^ for q"'" = -0. 135 and -1.35. The Nu^ for both cooling fluxes are simi- 
lar to results of constant property analysis for Gz“^ < 0.01. However, downstream of 
that point the results for the higher cooling flux diverge from those of the lower flux and 
never reach a constant value that can be characterized as "fully developed" heat transfer, 
although figures 4(a) and (b) indicate for q'*’ = -1. 35 a somewhat sinusoidal variation 
in NUj^ around the value 7. 0 ± 0. 4 in the larger Gz"^ range. 

On the other hand, the calculated Nuj^ for the q"'' = -0.135 develop closely along the 
curve for constant property fluids with initially uniform conditions (ref. 2). This classi- 
cal case asymptotically approaches the analytical value for the fully developed uniform 
heat flux case 


lim Nu = 

Gz“^— CO 


1 ®. 

11 


(ref. 5, p. 375). 

The Nujj increase for the q"*" = -1. 35 value, as compared with the q"*" = -0. 135 
development, is attributed to the steeper temperature wall gradient. Recall the second 
expression for Nusselt number, equation (29), 






or 
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The wall gradient (90/9R)j^_j^ increases quite rapidly in the region 0. 01 < Gz“^ < 0. 10, 
while the temperature decreases at a lower rate until Gz"^ = 0. 10, where the 

absolute wall temperature has reached a value low enough to cause the temperature ratio 
to counteract the increasing temperature gradient. At this point the Nuj^ reaches a 
maximum value of about 7. 0 to 7.2. The decrease beyond the maximum Nusselt is 
caused by very low calculated wall temperatures. The numerical program was halted in 
this region of reciprocal Graetz numbers. 

In reference 1, the influence of rate of heat flux on the fully developed Nusselt num- 
ber was demonstrated by a correlation between Nu^^ and the ratio By means 

of cross plotting data of the present numerical analysis, the developing Nusselt numbers 
were obtained as functions of the local ratio. In the limit of fully developed heat 

transfer, the numerical Nuj^ will approach the values of reference 1 as a limit. Fig- 
ure 5(a) demonstrates that this does happen for the low heat flux case. Note that the 
data in figure 5(a) is divided into the entrance region of developing velocity and tempera- 
ture profiles and the downstream region where velocity and temperatures are almost 
fully developed. Since there is no sharp separation point between these two regions, 
this division in figure 5(a) has arbitrarily been placed where our numerical Nusselt value 
is 110 percent of the fully developed Nusselt number of 48/11 for constant property anal- 
ysis. 

Two numerical runs for q'*' = -1. 35, shown in figure 5(b) again demonstrate the 
marked effect of strong cooling on the boundary layer development and that effect, in 
turn, on derived variables such as Nusselt number. The reversal of the Nusselt number 
curve produces maximum values about 75 percent higher than the predictions from con- 
stant property analysis. Now up to Z'*’ ~ 0. 01 the similarity in development of velocity 
and temperature profiles for both q"^ = -0. 135 and -1. 35 (see figs. 2 and 3) is also 
observed in the developing Nusselt numbers (fig. 4). The developing flow region in 
figure 5(b) corresponds approximately to the same Z'*’ as for that region in figure 5(a). 
There is, however, no quasi -developed flow region in figure 5(b) at the lower values of 
the temperature ratio there one finds instead the effects of rapidly increasing 

flow distortion due to the high cooling flux. 


Friction Parameter and Pressure Drop 

Friction factors based on wall shear stress are principally a function of the velocity 
profile gradient at the wall. The friction factor is defined by equation (34), and the 
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modified friction parameter 4f^ equation (36). 

Figures 6(a) and (b) give friction parameter results from the present analysis com- 
pared with Langhaar’s linearized entrance region analysis (ref. 7) and the numerical 
boundary layer analysis for constant property fluids of reference 2. Figure 6(a) shows 
a very close similarity of development of the friction parameter between the present 
analyses for q“*“ = -0. 135 and the constant property calculations of reference 2. 

With wall cooling fluxes of = -1. 35, the velocity profiles with inflection points 
near the wall (figs. 2(c) to (e)) produce smaller wall shear stresses than for developing 
entrance flow of incompressible fluids. This in turn produces the lower friction param- 
eters for the present numerical analysis which, in figure 6(b), are also compared with 
the calculations of references 2 and 7. 

Reference 1 correlates the shear -stress parameter 6 with the temperature ratio 
tw/tb to demonstrate the influence of heat flux rate on wall shear stress. As with the 
preceding discussion on the Nusselt number, the present numerical data for developing 
shear parameters was replotted to obtain the developing 5 parameter correlated with 
tw/tb* The results compared with the fully developed predictions of reference 1 are 
presented in figure 7. 

Both figures 7(a) and (b) have the similar entrance region development up to about 
^ 0. 01 where velocity and temperature profiles approach fully developed values 
(figs. 2 and 3). This has been noted also in figures 5(a) and (b) for Nusselt numbers. 
Beyond Z"’' 0. 01 for the moderate q"^ = -0. 135 there is quasi-developed flow, but 

the two almost identical curves for initial Mach numbers of 0. 004 and 0. 06 approach a 
limiting value of 6 below the prediction of the fully developed analysis in reference 1. 
This divergence of shear results between boundary layer calculations and fully developed 
analysis, even for essentially incompressible flow, was first noted in reference 2. Heat 
flux, whether heating or cooling, and even at moderate values, has a more marked effect 
on friction characteristics of flow than it does on the heat transfer parameters. 

At the higher cooling flux of q“^ = -1.35, the developing shear parameters 6 plotted 
in figure 7(b) are similar to those in figure 7(a) to around Z"*" ^ 0. 01. Downstream of 
that dimensionless length, the flow distortion due to cooling, which is discussed with 
reference to figures 2, 3, and 5(b), causes the shear stress to decrease rather rapidly 
with the lower values of the temperature ratio 

The behavior of the static pressure drop for gas flow through the tube is always of 
practical interest, and the pressure drop for flow with uniform cooling at the walls is 
of additional interest because of the adverse effect on the flow due to the energy with- 
drawal. Figures 8 to 10 describe the effects of the cooling on the pressure or pressure 
related functions. 

Figure 8 shows the average friction parameter (4fpbf^%)> based on the local static 
pressure drop from the entrance, as a function of the dimensionless distance Z"*". The 
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friction factor is given by equation (32), and the friction parameter actually plotted in 
figure 8 is given by equation (33). The friction factor (or parameter) is an averaged 
factor over the length Z"^. 

Figure 8(a) compares the present numerical results for moderate cooling and 
M. = 0. 004 and 0. 06 with results of constant property fluids from boundary layer anal- 
ysis (ref. 2), a linearized entrance region analysis (ref. 7), and the classical fully 
developed Poiseuille case. The present results are essentially the same as the earlier 
numerical results. 

The friction parameters for the higher cooling rate case at = 0. 004, 0. 06, and 
0. 30 are plotted on figure 8(b), Comparison is again made with the adiabatic incompres- 
sible numerical results of reference 2 and the linearized entrance analysis of reference 7. 
In addition, friction parameters shown in figure 8(b) extending from Z"*'= 0. 0001, have a 
curve closely fitting the experimental data of reference 8, which was for adiabatic incom- 
pressible laminar boundary layers in the entrance region of a tube up to Z'*'= 0. 001. The 
strong downward turn of the curves for the present numerical results starts around Z'*' = 
0. 01, which coincides with the similar behavior of the curves in figure 7(b) for wall shear 
stress and indicates the relation between the two variables. In the present calculations 
for higher cooling flux, the static pressure did rise above the initial p., which resulted 
in negative friction parameters. Since these could not be shown on the logarithmic plots 
of figure 8, this behavior will be discussed on the next figure. 

The development of dimensionless static and frictional pressure drops are recorded 
in figure 9, The dimensionless static pressure P'*’ was calculated directly in the numer- 
ical program. The pressure drop from momentum change in the flow, is related to the 
bulk velocity by equation (50) and the pressure drop due to frictional effects is the dif- 
ference between these two, as is given by equation (48). 

In figures 9(a) and (b) the pressure plots for = 0. 004 and 0. 06 at the moderate 
q'*' = -0. 135 indicate a small positive momentum pressure due to the developing velocity 
profile in the entrance region of the tube. The cooling soon reverses the sign on the 
momentum pressure so that the static pressure is less than the friction pressure for the 
rest of the tube length. The momentum pressure is not a significant part of the overall 
static pressure in the moderate cooling cases, and it is seen that the friction pressure 
and static pressure are virtually the same over most of the tube length. 

The pressure drops in the higher cooling level are shown in figures 9(c) to (e). The 
static pressure is seen to reach a maximum in the vicinity of Z"*" = 0. 01 and decreases 
steadily from that region. The difference between the friction pressure and the static 
pressure shown is the negative momentum pressure. It is interesting that the friction 
pressure for Z'*' > 0. 015 increases very slowly compared with its earlier behavior. 

Since the pressure values are cumulative along the tube, this almost constant value of 
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friction pressure indicates a sharp decrease in wall friction, which we have already 
observed in figure 7. 

These variations in the dimensionless pressure P'*’ indicate a rise in static pres- 
sure and a thickening of the wall boundary layer in the downstream tube sections at the 
higher cooling flux. This adverse pressure gradient recalls a similar situation in ex- 
ternal flows where thickening of the boundary layer precedes separation. In the present 
case of internal flows, the static pressure rise is accompanied by a continual decrease 
in the bulk velocity and persistance of flow in face of static pressure increase is due to 
the fluid inertia. 


Rayleigh Flow Calculations 

The shear stress friction factors of figures 6 and 7 show that the standard incom- 
pressible or fully developed correlations do not describe the frictional behavior of highly 
cooled gas flows, but they do describe the development of the moderate cooling cases. 

We can inquire whether the momentum ’’friction'' or pressure can be predicted by use 
of the Rayleigh calculations for one -dimensional frictionless heat transfer. These 
Rayleigh calculations have been carried out and are compared in figure 10 with the mo- 
mentum pressure drop calculated from the present numerical boundary layer analysis. 

The results of the momentum calculations for the moderate cooling levels for 
M. = 0. 004 and 0. 06 are given in figure 10(a). Rayleigh momentum pressures with 
cooling are always negative from the tube entrance. The boundary layer momentum 
pressure, even with cooling, is slightly positive close to the entrance because of the 
establishment of the velocity profile. The general shapes of the curve pairs are similar. 

Logarithmic plots of Rayleigh flow - boundary layer flow momentum pressures for 
the higher cooling fluxes are given in figures 10(b) to (d). The nature of the logarithmic 
scale accentuates the discrepancy between the two results at the lower values of the Z'*'. 
In terms of the absolute values, the Rayleigh and boundary layer results do not differ 
greatly. At the limit of the Z"*" range the difference is 10 to 20 percent. 

The interesting point with these calculations at q '" = -1. 35 is that the Rayleigh flow 
results underpredict the momentum pressure for the larger ranges of Z"*^. This is again 
due to the change in the velocity profiles in these Z"*" ranges. Profiles that are deviat- 
ing from fully developed contours have a positive momentum pressure component which 
opposes the negative component arising from the bulk velocity decrease. 
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SUMMARY OF RESULTS 


The compressible momentum and energy boundary layer equations are numerically 
solved for subsonic flow of a monatomic gas in a circular tube with initial uniform 
velocity and temperature profiles (at Z"*" = o), and with uniform wall cooling at flux 
levels q'*' of -0. 135 and -1. 35. The results of the calculations are as follows: 

1. At the lower cooling flux, the axial development of velocity and temperature pro- 
files, friction factors and Nusselt numbers, closely follow classical results for constant 
property fluids. 

2. For the higher cooling flux cases, the more important results are: 

a. In the downstream tube section, the velocity profiles thicken and the flow is 
similar to that of unconfined flows approaching a separation point. 

b. The increasing transverse velocity accompanying the thickening boundary 
layer augments thermal convection and flattens the radial temperature profiles. 

c. The bulk Nusselt number rises to a value 75 percent greater than the fully 
developed value of 48/11 for constant property fluids. 

d. Decreasing wall shear depresses the local friction factor values to less than 
one -half the Poiseuille limit. 

e. Static pressure increases above ambient as the fluid slows down under cool- 
ing, and adverse pressure gradients occur in conjunction with the thickening boundary 
layer. 

f. Rayleigh flow analysis gives momentum pressure drops within 25 percent of 
those obtained from the numerical boundary layer solution. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 16, 1970, 

129-01. 
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APPENDIX - SYMBOLS 


a 

s 

D 

Vb 

^T, b 

G 

Gz 

h 

K 

k 

M 

Nu^ 

P"^ 

Pe 


^FR 

p+ 

Pr, 


P 

T 

P 

Q 

q, 

q 

R 
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w 

+ 


local sound velocity, Vy^t 
heat capacity at constant pressure 
diameter of tube 

friction factor based on pressure drop with density evaluated at bulk tempera- 
ture, (+1/2) (p. - p)y/j(z/D)PjjU^J 

friction factor based on local wall shear stress with properties evaluated at bulk 
temperature, ( 1/ 2)PjjU^ 

mass flux 

Graetz number, (z/D) 

local heat transfer coefficient, 

dimensionless thermal conductivity, 

gas thermal conductivity 

local Mach number, u/a 

local bulk Nusselt number, hD/k^ 

dimensionless pressure drop, (p^ - p)^pm? 

Peclet number, Pr^Rej^ 

dimensionless pressure drop due to friction, eq. (49) 
dimensionless pressure drop due to momentum, eq. (50) 
initial Prandtl number, C u. /k. 

hr 

static pressure 

local static pressure ratio, p/p^^ 
cooling rate, energy per unit time 

cooling flux at wall, energy per unit time per unit wall area 
dimensionless cooling flux parameter, 
dimensionless radius, 


Rer. 


Re- 


gas constant 

bulk Reynolds number, 

initial Reynolds number , Duj^Pj^ Pj 
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r radial distance 

tube radius 

T dimensionless total temperature, - 1/2 [jy - l)/q^M?U^ 

dimensionless static temperature, (1 - 
t static temperature 

t’ local static temperature ratio, t/t^ 

U dimensionless axial velocity, u/u^ 
u axial velocity component 

V dimensionless radial velocity, 2(v/u^)Rej 

V radial velocity component 

W flow rate, mass per unit time 

Z dimensionless axial distance, (z/D)/Re^ 

Z*^ modified axial distance, (z/D)/Re|^ 

z distance along tube from entrance 

y ratio of specific heats, C /C 

^ V 

6^ shear stress parameter, based on bulk properties, r^r^/ 

6 static temperature ratio, (t - - t^) 

II gas viscosity 

ii' dimensionless viscosity, [i/ fi^ 

p gas density 

local wall shear stress 

7^ average wall shear stress over tube length z 
Subscripts: 

b bulk average at local axial position 

i initial (Z - 0) condition 

j,k mesh positions in difference technique 
o centerline (r = 0) conditions 

w tube wall conditions 
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Bulk Nusselt number, Nu 



Temperature ratio, 


(a) Initial Mach number, 0.004; dimensionless cooling flux, -0. 135 (b) Dimensionless cooling flux parameter, -1.35 {high cooling flux), 

(low cooling flux). 

Figure 5. - Comparison of developing local Nusselt numbers with values from analysis of fully developed flow with property variation. 
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Friction 


ZJ 



Present analysis 

Ref. 7 

Ref. 2 



(a) Initial Mach number, 0.004 and 0.06; dimensionless cooling flux parameter, -0.135 
(low cooling flux). 



Figure 6. - Friction parameter based on local shear stress plotted against dimension- 
less distance from tube entrance. 
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Shear stress parameter, 




(b) Dimensionless cooling flux parameter, -1.35 (high cooling flux. 


Figure 7. - Comparison of developing local shear stress parameter 
with values from analysis of fully developed flow with property 
variation. 
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Friction parameter, 4fp t|Re |3 = 




(b) Dimensionless cooling flux parameter, -1. 35 (high cooling flux). 

Figure 8, - Average friction factor based on static pressure drop plotted against 
dimensionless distance from tube entrance. 
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Friction pressure jj 


-Static pressure 


I 1 1 1 1 1 J I I 1 1 1 


(a) Initial Mach number, 0.004; dimensionless cooling flux parameter, -0.135 
(low cooling flux). 


Friction pressure--^ // 

y 

/A- static 


I L 1 1 iL 

io'2 10'^ 

Dimensionless axial distance, Z'^ 

(b) Initial Mach number, 0.06; dimensionless cooling flux parameter, -0.135 
(low cooling flux). 

Figure 9. - Static and frictional pressure drop in entrance length of tube 
with cooling. 


Friction pressure 



(c) Initial Mach number, 0.004; dimensionless cooling flux parameter, -1.35 (high cooling flux). 



(d) Initial Mach number, 0.06; dimensionless cooling flux parameter -1.35 (high cooling flux). 



Dimensionless axial distance, Z‘*‘ 

(e) Initial Mach number, 0.3; dimensionless cooling flux parameter, -1.35 (high cooling flux). 


Figure 9. - Concluded. 



Dimensionless momentum pressure drop, 



(a) Dimensionless cooling flux parameter. -0. 135. 

Figure 10. - Momentum pressure change by boundary layer calculation and Rayleigh flow (cooling). 
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\\ 

Rayleigh 

flow \ 


- Boundary layer 


lO'"^ 10'^ 10^ 

(b) Initial Mach number, 0.004; dimensionless cooling flux parameter, -1.35 (high 
cooling flux). 


\ \ Boundary layer 

■ \ \ 

\ \ 

Rayleigh, 

flow-'^ \\ 

\\ 

\ 


(c) Initial Mach number, 0.06; dimensionless cooling flux parameter, -1.35 (high 
cooling flux). 


- \ 

_ \ 

^ \ 

- \ 
“ \ 
Rayleigh 


-Boundary layer 


: I I I I I I II! ^ II I I n il l__i 1 I 1 1 1 1 r I I 1 I I 1 1 1 

^ 10 '^ 10 10 ’^ 1 
Dimensionless distance, Z'*’ 

(d) Initial Mach number, 0.3; dimensionless cooling flux parameter, -1,35 (high 
cooling flux). 

Figure 10. - Concluded. 
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